@myloaddata.pro
dxx=bdxx+2.0*bdxy+bdyy
;dxx=alog10(bdxx+2.0*bdxy+bdyy)
nstep=50
step=(max(dxx)-min(dxx))/nstep
level=min(dxx)+findgen(nstep)*step
xr= min(alpha)+findgen(4)*(max(alpha)-min(alpha))/3.0
yr= min(ek)*(max(ek)/min(ek))^((findgen(4)+1)/3.0)
;zr= alog10(min(dxx)*(max(dxx)/min(dxx))^((findgen(4)+1)/3.0))
device, decomposed=0
;loadct, 33,ncolors=nstep+1
;tvlct, 255,255,255,0
ct=colortable(33)
c=contour(dxx,alpha,ek,/fill,n_levels=nstep,rgb_table=ct,$
  background_color='white',/ylog,xstyle=1,xtickvalues=xr,$
  position=[0.1,0.1,0.85,0.95])
;contour,dxx,alpha,ek,/fill,nlevels=nstep,c_colors=reverse(indgen(nstep)+1),$
;  xstyle=1,ystyle=1,zstyle=1,/ylog
;cb=colorbar(target=c)
cb=colorbar(target=c,orientation=1,$
  textpos=1,position=[0.87,0.1,0.9,0.85],taper=0,$
  tickvalues=0,hide=0)
c.Save, '3.eps'
end